Astronomy & Astrophysics manuscript no. paper 


© ESO 2008 


February 5, 2008 





New evidence on the origin of the microquasar GRO J1 655-40 

J. A. Combi 1 , L. Bronfman 2 , 1.F. Mirabel 3 ' 4 

1 Departamento de Fi'sica (EPS), Universidad de Jaen, Campus Las Lagunillas s/n, 23071 Jaen, Spain 
e-mail: jcombi@ujaen.es 

2 Departamento de Astronomi'a, Universidad de Chile, Casilla 36-D, Santiago, Chile 
e-mail: leo@das.uchile.cl 

3 European Southern Observatory, Alonso de Cordova 3107, Santiago, Chile 

4 On leave from CEA-Saclay, France 
e-mail: fmirabel@eso . org 

Received; accepted 

ABSTRACT 

Aims. Motivated by the new determination of the distance to the microquasar GRO J1655-40 by Foellmi et al. (2006), we conduct a 
detailed study of the distribution of the atomic (HI) gas, molecular (CO) gas, and dust around the open cluster NGC 6242, the possible 
birth place of the microquasar (Mirabel et al. 2002). The proximity and relative height of the cluster on the galactic disk provides a 
unique opportunity to study SNR evolution and its possible physical link with microquasar formation. 

Methods. We search in the interstellar atomic and molecular gas around NGC 6242 for traces that may have been left from a supernova 
explosion associated to the formation of the black hole in GRO J1655-40. Furthermore, the 60/100 yum IR color is used as a tracer of 
shocked-heated dust. 

Results. At the kinematical distance of the cluster the observations have revealed the existence of a HI hole of l°.5xl°.5 in diameter 
and compressed CO material acumulated along the south-eastern internal border of the HI cavity. In this same area, we found extended 
infrared emission with characteristics of shocked-heated dust. Based on the morphology and physical characteristics of the HI, CO 
and FIR emissions, we suggest that the cavity in the ISM was produced by a supernova explosion occured within NGC 6242. The 
lower limit to the kinematic energy transferred by the supernova shock to the surrounding interstellar medium is ~ 10 49 erg and the 
atomic and molecular mass displaced to form the cavity of ~ 16.500 M G . The lower limit to the time elapsed since the SN explosion 
is ~ 2.2xl0 5 yr, which is consistent with the time required by GRO J1655-40 to move from the cluster up to its present position. The 
observations suggest that GRO 11655-40 could have been born inside NGC 6242, being one of the nearest microquasars known so 
far. 
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1 . Introduction lar atmosphere. The relative abundances of these elements sug- 
gest that the supernova progenitor could be a massive star with 

Microquasars are X-ray binary systems exhibiting relativistic ra- 25 to 40 solar masses 

dio jets (Mirabel & Rodriguez Egg. These systems contain Usi radio continuum and ffl -line observations around the 

compact objects like stellar black holes or neutron stars that ac- iousl a ted distance of 3 2 k for GRO J165 5_4o 

crete matter from a mass donor companion star. The formation (wWch corresponds t0 a velocity range from _ 42 t0 _ 30 km s -i )f 

of stellar black holes may result from direct collapse of mas- | C ombietal. (20011 found a cavity with "horse-shoe" morphol- 

sive star progenitors or from delayed collapse ofproto-neutron 'o g y in the large-scale HI distribution. This kind of structure has 

stars that undergo energetic explosions (Mirabel [2004|t. In some been already observed towards ea rly-type stars with very strong 

cases, compact objects can suffer a strong natal kick imparted in stellar winds (Be ha & c 1999) and supernova remna nts 

the supernova event. Therefore, they can escape from the region (g utt et a j 2001 ) 

of birth, where there could still exist signatures of their forma- u . , u . ,^r>/~» T i <rcc Ar, a 

" Other astrophysical connection between GRO J1655-40 and 

its environment has been suggested by Mirab el et al. (2 002). 
Until now, only 15 microquasars are known in the galaxy These authors haye used Hubble Space Telescope measuremen ts 

(Paredes 2005). Among them, GRO J1655-40 is one of the t0 compute a runaway space velocity of 1 12+18 km s" 1 for the 

most studied ones, with observations available along the en- object5 which is moving in opposite dire ction of the open cluster 

tire electromagnetic spectrum. The object is a binary system lo- NGC 6242 , being the latter located at {U b) = (345047, T A5). 

cated at (/, b) = (344!98, 2M5), which contains a 5-7 M©black According t0 photometric and radial velocity measurements 

hole with a companion star of 2.3 M© (Bailyn et al. [1995l>. Its (G l U shkova et al. I2006T> this open cluster is located at a dis- 

ongin has been disc ussed in the past in different frameworks. tance of L02±0 .l kpc from the Sun. Therefore, such connec- 

|Israelianetal.( 1999)| detected large overabundances of O, Mg, tion could not be confirmed since the old estimate of the dis- 

S i and S in the atmosphere of the companion star, which were tance for the microquasar had a value d i ffere nt from the one of 

interpreted as evidence for supernova ejecta captured by the stel- the open cluster Recently, the distance to GRO J1655-40 has 

been revisited by Foellmi ^t al. (2006*)| Using ESO archive and 

Send offprint requests to: J.A. Combi new VLT-UVES spectra they have determined a spectral type 
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F6IV for the secondary star of the binary system. Thus, the dis- 
tance to GRO J1655-40 would be smaller than 1.7 kpc. More 
recently, Foellmi (2006) showed that although this upper limit 
is rather firm, the lower limit is not well-constrained. Therefore, 
GRO J 165 5 -40 could be one of the closest known black holes 
to the Sun. 

In a recent paper on GRO J1655-40 by |Sala et al. (2007)| us- 
ing XMM data the authors derive a link between the inner radius 
of the accretion disk and the distance. They argue that if the dis- 
tance were smaller than 1.7 kpc and the mass is 5 Mq, then the 
inner accretion disk radius will be inside the gravitational radius 
of the black hole. However, the mass and spin of the black hole 
are uncertain. If the black hole were rapidly rotating as suggested 
bv lMcClintock et al. 20061 and its mass less than 5 Mq, the hori- 
zon of the Kerr black hole could be well inside the inner accre- 
tion radius. Furthermore, the authors admite that an estimate of 
the inner accretion radius using XMM data is model dependent. 

Motivated by the new determination of the distance to the 
microquasar, we conduct here a detailed study of the distribu- 
tion of the atomic (HI) and molecular (CO) gas, as well as of the 
dust in the vicinity of the open cluster NGC 6242. Our main goal 
is to detect the traces in the interstellar medium (ISM) around the 
cluster (i.e., a cavity), left by some explosive event. The proxim- 
ity and relative height of the cluster on the galactic disk provides 
a unique opportunity to study SNR evolution and the possible 
physical link with microquasar formation. The paper is orga- 
nized as follows: in Sect. [2 we describe our HI and CO observa- 
tions. In Sect. [3] we present the main results and then we discuss 
the nature of the detected structures, and their possible physical 
connections with GRO J1655-40, in Sect.ffl 

2. HI and CO observations 

2.1. HI observations 

We have performed HI observations towards GRO J 165 5 -40 
with a 30-m single dish telescope located at the Instituto 
Argentino de Radioastronomia (IAR), Villa Elisa, Argentina. 
The observations were carried out during six consecutive ses- 
sions on November 20-25, 2004. The radiotelescope has a 
helium-cooled HEMT amplifier with a 1008-channel autocor- 
relator, attaining a system temperature r sys « 35 K. System pa- 
rameters and a dditional details of t he observational technique 
can be found in Com bi et al. (1998)] The HI line was observed 
in hybrid total power mode and the sky was sampled on a 0.35° 
rectangular grid. Each grid position was observed during 60 s 
with a velocity resolution of ~ 1 km s" 1 and a coverage of ±450 
km s" 1 . A set of HI brightness temperature maps (Ar rms ~ 0.2 
K) were made for the velocity interval Av =(-80 km s" 1 , +20 
km s" 1 ). The data were analyzed using AIPS routine programs. 

2.2. CO observations 

The 12 CO data were observed during 2002, October with the 4 
m NANTEN telescope of Nogoya University, located then at Las 
Campanas (Carnegie Institution of Washington) in Chile. At the 
12 CO (7=1-0) transition of 115.27 GHz, the angular resolution 
(HPBW) is 2'. 8. A field of 1°.45 x 1°.45 centered on /=345°.45, 
b=+2°.5 was observed in the 12 CO line with an integration time 
of 40 s for each position. The typical rms noise antenna temper- 
ature of the observations is in the range of 0.2-0.25 K, and the 
velocity resolution is of 0.055 km s" 1 . 

The data were reduced with the NANTEN data reduction 
software. Baselines were fitted to each individual spectrum. The 
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Fig. 1. HI brightness temperature channel maps obtained for the 
velocity range from -14 to -6 km s" 1 around the position of 
NGC 6242 (indicated with a white small circle). In all images 
the contours are labeled in steps of 4 K, starting from 4 K. A 
J2000 frame is shown for reference in the figure and all images. 
The location and motion direction of GRO J 165 5 -40 are shown 
with a small circle and an arrow, respectively. 

resulting data cubes were transformed into FITS format, and fur- 
ther processing was done using the AIPS software. 

3. Main Results 

If the physical link between GRO J1655-40 and NGC 6242 is 
real, the system should have been ejected from the cluster ~ 3 Ax 
10 5 years ago. After this time the radio emission of a SNR could 
be very faint or not detectable. In order to confirm this hypoth- 
esis, we search for signatures of some explosive event around 
NGC 6242 at radio frequencies using data from the 4.85-GHz 
(Condon et al. [T993b and 2.7-GHz surveys (Duncan et al. [T995l) , 
respectively. As expected, we found no traces of any SNR phe- 
nomena towards the region. 

3.1. Atomic component 

Standard galactic rotation models (Russeil 2003), corrected by 
exces s associated with part of the Carina arm (Alvarez et al. 
1990), indicate that for a kinematic distance in the range of 0.8 
to 1.6 kpc (the minimum distance to NGC 6242 and the max- 
imun distance to GRO J 165 5 -40), corresponds a velocity inter- 
val of -14 to -6 km s" 1 . Fig.l, displays HI brightness tempera- 
ture channel maps of the relevant interval. The development of a 
cavity centered at (/, Z?)=(345.67°,+2.68°) can be clearly seen in 
these maps. The angular size of this structure is ~ 1°.5 x 1°.5. 
Thus, at an average distance of 1.2 kpc, the linear size of the HI 
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Fig. 2. Integrated column density map obtained for the velocity 
range from -14 to -6 km s" 1 around the position of NGC 6242 
(indicated with a white circle). The contours are labeled in steps 
of 4 K, starting from 4 K. The location and motion direction 
of GRO J 165 5 -40 are shown with a small circle and an arrow, 
respectively. 



structure would be ~ 32 pc. The location of NGC 6242 and the 
HI cavitiy are indicated with white circles. 

The angular extension of the cavity increases from -14 to 
-10 km s" 1 . Between -12 to -10 km s" 1 , the hole reaches its 
maximun size and forms a kind of "horse- shoe" structure at ve- 
locities larger than -8 km s" 1 . The HI map integrated over the 
velocity interval from -14 to -6 km s" 1 is shown in Fig. 2. 

Given that, no radio emission was detected towards the HI 
hole, we have computed the upper-limit for the radio surface 
brightness emission at 1 GHz. Flux densities at 2.7 and 4.8 GHz 
of the region covered by the HI hole and an angular size for the 
cavity of 1.5° were used. As may be expected, it turns out to be 
very low, with £ 1G h z ~ 4(±2)xl0 -24 W m" 2 Hz" 1 sr" 1 . 

3.2. Molecular component 

Figure 3 shows the distribution of the CO(/=l— >0) emission in 
a smaller area around NGC 6242 (indicated with a black circle). 
A large amount of compressed CO lies along the lower half of 
the image, in the region closer to the Galactic plane. This ex- 
tended molecular feature coincides with the lower border of the 
HI cavity. It has a semi-circular morphology, typical of the inter- 
action of a SNR shock with the steep density gradient towards 
the galactic plane. From -16 to -6 km s" 1 the angular extension 
of the cavity increases with the velocity. It is evident that in this 
velocity interval the molecular material was displaced to form 
part of the cavity seen in the HI distribution. 




tion of NGC 6242 (indicated with a black circle). White dashed 
circles show the expansion of the CO material along the velocity 
range. The location and motion direction of GRO J 165 5 -40 are 
shown with a small circle and an arrow, respectively. 



At the south region of the cluster the CO distribution exhibits 
a peculiar bifurcation. The morphology of the material could be 
the result from the eruption of a shock into a region of very 
low density. When the shock front hits a low density region it 



Combi et al.: New evidence on the origin of the MQ GRO J 1655 -40 



04 oo 



03 30 



S 

u 

1 
i 



02 30 




03 10 



01 30 



347 00 



346 30 



00 345 30 

GALACTIC LONG. 



344 30 



Fig. 4. Comparison of the HI distribution (in grey- scale) in the 
velocity range from -14 to -6 km s" 1 with the 12 CO emission 
integrated in the same velocity range. In the HI image the con- 
tour are labeled in steps of 4 K, starting from 4 K. The plot- 
ted 12 CO contours (white lines) are 2, 2.5, 3, 3.5, 6.5, 9.5, 16.5, 
25.5 and 35.5 K km s" 1 . The location and motion direction of 
GRO J 165 5 -40 are shown with a small circle and an arrow, re- 
spectively. 



expands more rapidly into the interstellar medium producing a 
separation. The good agreement of the HI material, similar mor- 
phology and behavior of the molecular gas in this velocity range, 
suggests that all these structures have a physical link and origi- 
nate in the same explosive event. 

A comparison of the HI distribution integrated in the velocity 
range from -14 to -6 km s" 1 with the 12 CO emission integrated 
in the same velocity range, is shown in Fig. 4. It is interesting 
to note that the molecular material is mainly displaced towards 
the lower region of the HI cavity. This is consequence of the 
fact that the mean thickness of the molecular disk is less than 
that of atomic gas (Bronfman et al. 119881) . The circle indicate 
the position and aproximate extension of the open cluster NGC 
6242. 



3.3. Infrared emission 

We have also searched for traces of a SNR shock using IRAS and 
MSX observations. Figure 5 shows the MSX infrared emission at 
8.28 yum around NGC 6242. The size of the image represents the 
HI hole shown in Fig. 2. In spite that is a very confuse region of 
the Galaxy, we found that several IR structures with extended 
morphology have good correlation with the molecular emission 
along the southern part of the image. At the central region of 
the HI hole we found extended IR emission with a semi-circular 
morphology (feature A). For all these structures the IR emission 
follows the spectral trend F12 < F25 < F60 < F100 (where F^ is the 
flux density at wavelenght A) characteristic of shocked-heated 
dust (Junkes et al. 1992 ). At the southern part of the image, we 
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Fig. 5. Infrared emission from MSX at 8.28 fim. The position of 
the NGC 6242 is indicated with a black circle. The location of 
the center of the HI hole is indicate with a black cross. The size 
of the figure is approximately that of the HI cavity. The central 
semi-circular feature is labeled with letter A. 



can also see a large amount of IR emission, which could repre- 
sent the shocked dust component of the associated CO material. 

An increase in the 60/100 jim IR color is a useful tracer of 
shocked-heated dust. The 60/100 jim color-corrected ratio for 
the observed infrared emission of the whole region is shown in 
Fig. 6. It is interesting to note, that the infrared ratio level changes 
abruptly on the south border of the HI hole, which suggests the 
presence of the SNR shock. The small black circle represents 
the aproximate extension of NGC 6242. Enhanced IR emission 
on the central region and the southern part of the HI hole, is 
revealed by the ratio 60/100 fim. This values agree very well 
with the color temperature derived for IRAS sources associated 
with shocked material in old SNRs (Huang et al. 1986). 

4. Discussion 

Atomic and molecular observations reveal that NGC 6242 is sur- 
rounded by shocked material with semi-circular structure, typi- 
cal of a medium where a large amount of energy has been re- 
lease. The existence of a local minimum in the HI distribution 
at the approximated position of the cluster is evident. It is con- 
firmed by the CO distribution which is bending on the lower bor- 
der of the HI cavity. Moreover, infrared emission found within 
and along the south-eastern border of the HI hole, has proper- 
ties typical of dust heated by collisions with the postshock gas. 
All these facts suggest that the detected emission has an under- 
lying physical connection originated in the same explosive phe- 
nomenon. 

The energetic impact of a SN explosion that occurs out of 
the galactic disk can drive an asymmetric expanding bubble if 
the surrounded ambient gas is non-uniform. The shock-front 
could severely affect the dynamic structure of the surrounding 
gas, originating the observed semi-circular cavity, and heating 
the ambient dust. Assuming that the HI - CO cavity is the signa- 
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Fig. 6. Gray-scale distribution of the 60/100 jim infrared ratio, 
which varies between 0.5 and 5. Note as the infrared ratio level 
changes abruptly on the south border of the HI hole. The position 
of the HI hole and NGC 6242 are indicated with white circles. 
The location and motion direction of GRO J 165 5 -40 are shown 
with a small circle and an arrow, respectively. 



ture left by the action of a SNR on the surrounding ISM, we can 
use the N H I images and CO maps to estimate the mass swept- 
up by the front- shock. It allows to obtain the initial ambient gas 
density and roughly to estimate the time elapsed since the SN 
explosion. 

Estimates of the physical parameters for the HI structure 
were derived by fitting a circle to the maximum of the envelope, 
which covers an area ofA/xAZ? « 1°.5 x 1°.5. The integrated 
column densities shown in Fig. 2 were used to estimate the to- 
tal mass displaced to form the cavity. If the background density 
is given by the density at the minimum, we obtain that at least 
9.000 M Q have been removed from this region. The expansion 
velocity of the shell is defined approximately as (Av/2)+2 km 
s" 1 , where A=V2-vi is the velocity interval through which the 
cavity is seen. Using an average expansion velocity of -12 km 
s" 1 , we obtain a lower limit to the original energy released in the 
explosion of ~ 2xl0 49 erg. Since the SNR can have lost a sig- 
nificant part of its energy via radiative cooling, the initial energy 
could have been ~ 10 51 erg. The mean density before the explo- 
sion (i.e. the initial density) can be obtained using the removed 
mass from the envelope's volume. Averaging the column densi- 
ties within the cavity, we obtain that the density due to atomic 
material is ~ 14 cm" 3 . 

We further estimate the molecular component (H2) of the 
initial density, averaging 41 molecular clouds of the Carina arm 
taken from the study of Grabelsky et al. (1988) Taking into ac- 
count the region of the lower half of the HI hole and including 
the correction of He, we obtain that the density of molecular 
hydrogen would be ~ 14 cm" 3 , that corresponds to 7500 M© 
(which renders a total atomic and molecular displaced mass of 
16.500 M Q ). Assuming that the whole cavity has been emptied 
of HI material and only the lower half of CO, we estimate that 



an upper limit for the initial density, averaged over the whole 
cavity, is N(H) ~ 28 cm" 3 . With this density and a radius of 16 
pc for the cavity, we can estimate the lower limit for the SNR's 
age using the dynamical evolution equations for an expanding 
SNR, described by |Koo & Kang (2004)[ The non-detection of 
radio emission and size of the SNR's radius, suggest that it is 
well within the radiative phase. The transition to this phase in 
SNR evolution is often triggered when the blast wave interacts 
with dense atomic and molecular clouds of gas in the ISM, a 
fact supported by the large amount of shocked gas observed in 
the region. 

According to [Koo & K ang (2004) , the maximum age of the 
HI hole due to the SNR is defined as the time when the shock ve- 
locity drops to the ambient sound speed, so that the interior pres- 
sure is comparable to the ambient pressure and the shell breaks 
up. In this case, the lifetime of the SNR is described by: 

r = 2.02 x 10 6 OSc s , 6 )" L43 n- a367 ^ 1 316 yr - 0.91t sf (1) 

where c^6=c^/(10 6 cm s" 1 ), c^ is the sound speed, /3=2, E51 
is the SN energy released to the ISM in units of 10 51 erg, no is 
the medium density, and t S f is the shell formation time given by, 

t sf = 3.61 x 10 4 rc~ 4/3 £;?( 14 yr. (2) 

Using a mean density of no -28 cm" 3 and a SN energy of ~ 
10 51 erg, we obtain t s f= 400 yr. For c Sr e=l, the lower limit to the 
lifetime of the SNR is r= 2.2xl0 5 yr. 

This value is on the order of the time elapsed by 
GRO J1655-40 in moving from NGC 6242 up to its present 
position, which is ~ 3.4x 10 5 yr. 

5. Summary 

We have investigated the surroundings of the open cluster NGC 
6242 through HI, CO, and IR observations. Several facts corrob- 
orate the hypothesis of a SN event occurring within or very close 
to NGC 6242. 

- The HI brightness temperature channel maps for the distance 
range of 1.2+0.4 kpc show the development of a cavity with 
an angular size of ~ 1°.5 x 1°.5. 

- The distribution of the CO(/= 1 -0) emission displays a large 
amount of compressed material with semi-circular morphol- 
ogy along the south and east borders of the HI cavity. 

- Infrared observations reveal the presence of several IR struc- 
tures with extended morphology and good correlation with 
atomic and molecular emissions. The 60/100 yum color- 
corrected ratio shows also enhanced infrared emission to- 
wards the south-eastern part of the HI hole. These values 
agree very well with the color temperature derived for IRAS 
sources associated with shocked material in old SNRs. 

- All the available information suggests that the detected emis- 
sions morphologies have a physical link, originated in a su- 
pernova explosion occurred within or very close to the clus- 
ter a tk 2.2x 10 5 years. If the association between the open 
cluster and the SNR is confirmed then the distance of the 
microquasar would be better known. 

Now, is the origin of GRO J 165 5 -40 physically associated to 
the explosive event?. We suggest that a possible scenario, which 
fits the data presented in this work, is that an old supernova ex- 
plosion occurred at least 2.2x 10 5 years ago inside NGC 6242. In 
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this explosive event GRO J1655-40 could originate and conse- 
quently be ejected up to its present position in the sky. A similar 
case, not confirmed, could be that of the m icroquasar LS 5039 
and the SNR G016.8-01.1 (IRibo et al. 2002b . 

Unveiling the origin of the microquasar will allow to con- 
strain its distance. If the present evidence in this work is con- 
firmed, GRO J 165 5 -40 would be the first microquasar physi- 
cally linked to a very dilute, but still detectable SNR. 

High-resolution HI observations are required for a most de- 
tailed study of the cold gas and evidence of cloud compression 
originated by the shock front. Investigations of other molecular 
species, such as OH, HCO+ and H2CO, could help in search of 
shock signatures, particularly around the south-eastern region of 
the HI hole, where a large amount of gas and dust is detected. 
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